Perdurance of multiply connected de Sitter space 
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This paper deals with a study of the effects that spherically symmetric first-order metric pertur- 
bations and vacuum quantum fluctuations have on the stability of the multiply connected de Sitter 
spacetime recently proposed by Gott and Li. It is the main conclusion of this study that although 
such a spacetime is stable to the classical metric perturbations for any size of the nonchronal region, 
it is only stable against the quantum fluctuations of vacuum if the size of the multiply connected 
region is of the order the Planck scale. Therefore, boundary conditions for the state of the universe 
based on the notion that the universe created itself in a regime where closed timelike curves were 
active and stable, still appear to be physically and philosophically so well supported as are those 
boundary conditions relying on the notion that the universe was created out of nothing. 

PACS number(s): 98.80.Cq, 04.20.Gz, 04.25. Nx 



I. INTRODUCTION 

When he was arguing in support of his favourite con- 
dition, Hawking considered [1] that the choice of the 
boundary conditions for the quantum state of the uni- 
verse should be done on just compact metrics or noncom- 
pact metrics which are asymptotic to metrics of maximal 
symmetry. While completedness seems to be an appeal- 
ing property, spaces which are maximally symmetric are 
wanted for their ellegance and greatest simplicity. At the 
time being, the quest for the boundary conditions of the 
universe has already collected a number of proposals with 
the above or similar properties, including the Vilenkin's 
tunneling wave function [2], the vanishing Weyl-tensor 
condition of Penrose [3] and, perhaps most popularly, the 
no boundary proposal of Hartle and Hawking [4] . They 
have all been much discussed from different perpectives 
[5], particularly by invoking their capability to generate 
a suitable inflationary mechanism able to solve original 
problems of standard cosmology [5,6], giving at the same 
time rise to a scale invariant spectrum of density fluc- 
tuations [7,8], compatible with the present observational 
status of the universe. 

Recent observations on the density of gravitational 
lenses in the universe [9,10] and estimates of the present 
value of the Hubble constant [11,12] have nevertheless 
led to the conclusion that the critical density £1 has to 
be smaller than unity, so favouring an open rather than 
closed model for the universal expansion. The problem 
then is that most of the inflationary models compatible 
with the proposed boundary conditions predict a value 
for Q very close to unity. In order to solve this problem, 
Hawking and Turok have recently suggested [13] the ex- 
istence of a singular instanton which is able to generate 
an open inflationay universe within the framework of the 
no boundary proposal. Vilenkin has argued [14] that the 
singularity of this instanton may have catastrophic con- 
sequences and, therefore, some procedures have been ad- 



vanced [15,16] to make the instanton nonsingular, while 
still inducing an open inflationary process. 

It is worth noting that so early as 1982 Gott already 
proposed [17] a procedure leading to an open inflation- 
ary model for the universe. At the time, of course, it was 
largely disregarded, since it was then generally believed 
that the universe is closed. Motivated by the present ob- 
servational status, Gott and Li have now suggested [18] 
new boundary conditions based on the Gott's original 
model of open inflation. These boundary conditions as- 
sume the existence of a nonchronal region with closed 
timelike curves (CTC's) in de Sitter space, separated 
from the observable universe by a chronology horizon 
[19]. The resulting model can be confronted to e.g. the 
most popular no boundary paradigm as follows. Con- 
sider first de Sitter space and visualize it by means of the 
Schrodinger five-hyperboloid with Minkowskian coordi- 
nates [20,21]. If we slice this space along surfaces of con- 
stant timelike coordinate, then the slices become three- 
spheres and represent a closed universe. On this slicing 
both the no boundary and the Vilenkin conditions hold, 
so that in these models the beginning of the universe can 
be pictured (in the Euclidean framework) as being the 
south pole of the Earth [22] , with the time running along 
the Earth's meridians. Thus, asking what happened be- 
fore the beginning of the universe is like asking what 
is south of the south pole. Therefore, these boundary 
conditions can be regarded to implement the idea that 
the universe was created out of nothing. In contraposi- 
tion to this approach, one can also slice the Schrodinger's 
five-hyperboloid vertically, i.e., along the only spacelike 
direction defined in terms of the proper time [18]. The 
resulting slices are negatively curved surfaces, describ- 
able in terms of open cosmological solutions. It is on this 
slicing that the Gott-Li's boundary condition can be de- 
fined. Since the time should now lie along the Earth's 
Ecuator (or actually any of its parallels), the origin of 
time can by no means be visualized or fixed, so asking 



1 



what was the earliest point is like asking what is the east- 
ernmost point on the Earth's Ecuator: there will always 
be an eastern (that is, earlier) point. This was implied to 
mean [18] that on some region of the slices there must be 
CTC's. These curves will make the job of always shifting 
the origin of the cosmological time. Saying then that the 
universe was created from nothing would be meaningless; 
what one should instead say is that it created itself [18]. 

We have thus two different types of boundary con- 
ditions of the universe that can induce it to be open. 
Whereas the no boundary condition does it by rather 
an indirect way which involves some suitably modified 
version of the Hawking- Turok instanton, the Gott-Li's 
proposal creates the open universe directly. Is the latter 
proposal therefore more fashionable than the former?. A 
positive answer to this question could only be made once 
the Gott-Li's model would satisfactorily pass some im- 
portant tests on its consistency. First of all, one had 
to check whether the multiply connected de Sitter space 
is classically and quantum-mechanically stable. Li and 
Gott claimed [23] that all multiply connected spacetimes 
with a chronology horizon (derived from Misner space) 
are stable to quantum fluctuations of vacuum, but pre- 
vious work by Kay, Radzikowski and Wald [24] and by 
Cassidy [25] has raised compelling doubts on this con- 
clusion. The present paper aims at partly filling the 
above requirements by studying the classical and quan- 
tum stability properties of the multiply connected de Si- 
iter space. This will be done using both, a first-order 
perturbation procedure paralleling the method devised 
by Regge and Wheeler [26] to investigate the stability of 
Schwarzschild spacetime, and a time-quantization proce- 
dure [27] to analyse regularity of the solution against vac- 
uum quantum fluctuations. Our main conclusion is that 
the multiply connected de Sitter space is stable both clas- 
sically and quantum-mechanically. Quantum stability is 
however restricted to hold only on the very small regions 
where the time shows its essential quantum character. 

We outline the paper as follows. In Sec. II we briefly 
review how a multiply connected de Sitter space can be 
constructed, and show why its nonchronal region must be 
confined inside the cosmological horizon. The stability of 
the whole de Sitter space has already been investigated 
in terms of a global Fricdmann-Robertson- Walker met- 
ric [28,29]. However, as far as I know, no corresponding 
research has been hitherto attempted for the de Sitter 
region covered by static coordinates, that is the region 
where multiply connectedness and CTC's should appear. 
We have performed this study here, first for the simply 
connected static case in Sec. Ill, and then for the mul- 
tiply connected case in Sec. IV. This Section also con- 
tains the analysis of the stability of multiply connected 
de Sitter space against quantum fluctuations of vacuum. 
Finally, we summarize and conclude in Sec. V. 



II. THE MULTIPLY CONNECTED DE SITTER 
SPACE 

de Sitter space is usually identified [21] as a maximally 
symmetric space of constant negative curvature (positive 
Ricci scalar) which is solution to the vacuum Einstein 
equations with a positive cosmological constant, A > 0. 
Following Schrodinger [20,29], it can be visualized as a 
five-hyperboloid defined by 

w 2 +x 2 +y 2 + z 2 -v 2 =p 2 , (2.1) 

where p a = ^/3/A. This hyperboloid is embedded in 
E 5 and the most general expression for the metric of 
the de Sitter space is then that which is induced in this 
embedding, i.e.: 

ds 2 = -dv 2 +dw 2 +dx 2 +dy 2 +dz 2 , (2.2) 

which has topology RxS A , invariance group SO (4, 1) and 
shows ten Killing vectors (four boosts and six rotations). 

Metric (2.2) can be coveniently exhibited in cither 
global or static coordinates. Global coordinates t' £ 
(—00,00), -03,-02 € (0, 7r), and 0i £ (0, 2ir) can be de- 
fined by [29] 

z = po cosh (t' I po) sin 3 sin 2 cos 0i 

y = po cosh (0 /po) sin -03 sin0 2 sin-0i 

x = p cosh (t'/po) sin-03 cos-0 2 (2.3) 

w — po cosh (0 / po) cos -03 

v = po sinh (t'/po) . 

In terms of these coordinates metric (2.2) becomes 

ds 2 = -dt' 2 + pi cosh 2 (t'/po) dnl, (2.4) 

where cftl§ is the metric on the unit three-sphere. Met- 
ric (2.4) is a k — +1 Fricdmann-Robertson- Walker met- 
ric whose spatial sections are three-spheres of radius 
po cosh(t' / po)- Coordinates (2.3) enterely cover the four- 
dimensional de Sitter space which would first contract 
until t 1 = and expand thereafter to infinity. 

In order to exhibit metric (2.2) in static coordinates 
t £ (—00,00), 03,02 € (0,7r), 0i £ (0, 2ir), one can use 
the definitions [29], 

z = po sin 03 sin 2 cos 0i 
y = Po sin 03 sin 2 sin 0i 
x = po sin 03 cos 02 (2.5) 
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w = po cos tp3 cosh (t/po) 

v = po cos ^ 3 sinh (t/p ) ■ 

Setting r = p sinip 3 (i.e. denning r e (0,po)) 5 we obtain 
the static metric in de Sitter space 

ds 2 = - (l - dt 2 (l - ^ dr 2 +r 2 dill (2.6) 

where dil 2 is the metric on the unit two-sphere. The 
coordinates defined by Eqs. (2.5) cover only the portion 
of the de Sitter space with w > and x 2 + y 2 + z 2 < p§; 
i.e., the region inside the particle and event horizons for 
an observer moving along r = 0. 

In order to see whether the whole or some restricted 
region of the de Sitter space can be made to have multi- 
ply connected topology, with CTC's on it, we will fol- 
low the procedure described by Gott and Li [18], so 
checking whether a symmetry like that is satisfied by the 
Minkowskian covering to Misner space [21] is somewhere 
holding in de Sitter space. On the Minkowskian five- 
hyperboloid visualizing de Sitter space, such a symmetry 
would be expressible by means of the identification [27] 

(v, w, x, y, z) (v cosh(n6) + w sinh(nfe), 

w cosh(nfe) + v sinh(nfe), x, y, z) , (2.7) 

where & is a dimensionlcss arbitrary quantity and n is 
any integer number. The boost transformation in the 
(v,w)-planc implied by this identification will induce a 
boost tranformation in the de Sitter space. Hence, since 
the boost group in de Sitter space is a subgroup of the de 
Sitter group, either the static or the global metric of de 
Sitter space can also be invariant under symmetry (2.7). 

It is easy to see that there cannot exist any symmetry 
associated with identification (2.7) on the (w,w)-palnc 
which leaves metric (2.4) invariant for coordinates (2.3). 
It follows that the whole of the de Sitter spacetime can 
neither be multiply connected, nor have CTC's. How- 
ever, for coordinates defined by Eqs. (2.5) leading to the 
static metric with an apparent horizon (2.6), the above 
symmetry can be satisfied on the region covered by such 
a metric, defined by w > \v\, where there are CTC's, with 
the boundaries at w = ±u and x 2 +y 2 + z 2 = p\ being the 
Cauchy horizons that limit the onset of the nonchronal 
region from the causal exterior [18]. Such boundaries be- 
come then appropriate chronology horizons for de Sitter 
space. 

III. STABILITY OF STATIC DE SITTER SPACE 

While our discussion of Sec. II made it clear that mul- 
tiply connected de Sitter space is mathematically rich 
and interesting, we still need to know if such a space 



is indeed a physical object. Therefore, in what follows 
we shall use a general-relativity perturbation method to 
investigate the stability of the multiply connected de Sit- 
ter universe. Since multiply-connectedness and CTC's 
only appear in the region covered by static coordinates, 
the extension of the analysis of the cosmologically per- 
turbed global metric for simply connected de Sitter space 
[28,29] to a multiply connected topology would unavoid- 
ably lead to rather inconclusive results. We instead shall 
proceed as follows. We first extend the perturbative pro- 
cedure originally devised by Regge and Wheeler for the 
Schwarzschild problem [26] to a cosmological de Sitter 
space in this Section, and then in the next Section, wc 
conveniently include the effects derived from the identi- 
fication (2.7) in the resulting formalism. In the present 
paper we confine ourselves to the linear analysis, investi- 
gating the stability of simply and multiply connected de 
Sitter space in first order perturbation theory by means 
of a generalization from the refined method developed by 
Vishveshwara [30] and Zerilli [31]. 

We take as the general background metric g M „ and the 
perturbation on it as h^. The quantity g^ will be later 
specialized to be the static de Sitter metric, i.e., 

ds 2 = - (l - H 2 r 2 ) dt 2 + (l - B 2 r 2 y X dr 2 + r 2 dVL 2 2 , 

(3.1) 

where we have now denoted H — pQ 1 for the sake of sim- 
plicity, and x° — t, x 1 — r, x 2 — 9, x 3 — <fr. Metric (3.1) 
corresponds to the initial time-independent equilibrium 
configuration, so the problem to be solved is, then, if 
metric (3.1) is somehow perturbed, whether the pertur- 
bations will undergo oscillations about the equilibrium 
state, or will grow exponentially with time. The static 
de Sitter space will be stable in the first case and instable 
in the second one. 

Since the background is spherically symmetric, any ar- 
bitrary perturbation can be decomposed in normal modes 
given by [26] 

^/o(f)/i(r)/ 2 (0)/ 3 W. 

Associated with these modes we have an angular moment 
I and its projection on the z axis, M. For any given 
value of I there will be two independent classes of per- 
turbations which are respectively characterized by their 
parities (— l) e (even parity) and (— (odd parity). 

Furthermore, since the background is time-independent, 
all time dependence of the perturbations will be given by 
the simple factor exp(— ikt), where k is the frequency of 
the given mode. 

In order to derive the equations governing the pertur- 
bations, we shall start with the Einstein equations 

R^ijg) = Kg^, (3.2) 

with g denoting the de Sitter background metric. For 
the perturbed spacetime, these field equations transform 
into 
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R^{g) + SRf, u (h) = A (g^ + h IMV ), 



(3.3) 



for small perturbation, with 8R^ V assumed to contain 
first-order terms in h^ v only. Now, since Einstein equa- 
tions are still valid in the perturbation scheme, we ob- 
tain in this approximation that the differential equations 
that govern the perturbations should be derived from the 
equations 



SR^(h) = Ah^ 



(3.4) 



The SR^v's will be here computed using the same general 
formulas as those employed by Regge and Wheeler [26] 
and Eisenhart [32]; i.e. 



(3.5) 



where the semicolon denotes covariant differentiation, 
and the variation of the Christoffel symbols is given by 

*T2„ = \g pa ( W + *W - V;«) • ( 3 - 6 ) 

After introducing suitable gauge transformations 
[30,31], the most general perturbations in static de Sitter 
space can be written in forms which are similar to those 
obtained for Schwarzschild space; i.e. 

For odd parity: 



h (r) 

hi(r) 

"0 

sym sym 



sinf 



no 



Pi(cos8) exp(-ikt). 



(3.7) 



For even parity: 



h 



H (1 - H 2 r 2 ) H, 

H, ^(l-i/V)" 1 







Kr 2 

Kr 2 sin 2 9 



horizon. In Kruskal coordinates u, v, the static de Sitter 
metric can thus be written as 

ds 2 = F(r) 2 (du 2 -dv 2 )+r 2 dtl 2 2l 

where 



(3.9) 



F{rf 



2 2 

V — U = 



1 + Hr 
H 



1 - Hr 
1 + Hr 



and 



Ht = - tanh 1 (- 



(3.10) 



(3.11) 



(3.12) 



In terms of the Kruskal coordinates, the components of 
the metric perturbations take the form (the components 
that invove only angular coordinates are the same in the 
two coordinate systems): 



h K 



F 2 



u 2 (l- H 2 r 2 ) h m 



+v 2 (1 - H 2 r 2 ) /in - 2uvh m ] 



h K - 



v 2 (1 - H'r 2 )- 1 h 00 



+u 2 (1 - H 2 r 2 ) /in - 2uuft i] 



l 01 — ^5 2 



[(u 2 +v 2 ) /ioi 



-uv 



((l-H 2 r 2 ) 1 /i 00 + (l-ff 2 r 2 )/iii)] (3.13) 



^03 a — 2 2 



[uh 03 -v(l- H 2 r 2 ) h 13 



x P £ (cos0)exp(-ifci). (3.8) 

In these expressions "sym" indicates that = h vtl , 
Pe(cos9) is the Legendre polynomial, and ho, hi, H , 
Hi, H2 and K are given functions of the radial coor- 
dinate r which must be determined as solutions to the 
respective wave equations, subject to suitable boundary 
conditions. In order to ensure holding of regularity on 
the cosmological horizon at r = H^ 1 , we should trans- 
form the components of the perturbations in Eqs. (3.7) 
and (3.8) to the representation where the static metric is 
maximally extended and regular also on the cosmological 



h i3 « 2 2 

u z — v 



[vh 03 - u (1 - H 2 r 2 ) /113] , 



where the superscript K refers to Kruskal coordinates. 
For future reference, we introduce here the relation: 



exp (Hr*) = u 2 — v 2 , 
where the new variable r* is defined by 



r*=H-^J Hr - 1 



\Hr + l 



(3.14) 



(3.15) 
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A. The wave equations for the perturbations 

We derive in what follows the differential equations 
which should be satisfied by the perturbations on the 
maximally extended de Sitter metric. We shall start with 
the odd-parity solutions, and choose the perturbed Ein- 
stein equations 



kh 



1-H 2 r 2 dr 



(f - H 2 r 2 ) hi=0 



(3.16) 



k % h HW ] + {e - m + 2) ^ - 247rGif2 ^ 

(3.17) 

where the prime denotes differentiation with respect to 
radial coordinate, = d/dr. Let us first introduce the 
definition 



(3.8), we obtain first the independent first-order pertur- 
bations of the Ricci tensor 

dK K-h H 2 rK ! i(£+l)h 1 _ 

dr + r + (1 - H 2 r 2 ) 2kr 2 ~ U [6 b) 



ikh x + (1 - H 2 r 2 ) - - 2H 2 rh = (3.26) 

v dr dr 1 



dr 



[(l-H 2 r 2 )h 1 ]+ik(h + K)=0 (3.27) 



2 d 2 h 2 



dh 



dr 



rQ = {l - H 2 r 2 ) hi, 
with which Eqs. (3.16) and (3.17) become 

— (rQ) =0 



[l-H 2 r 2 ) dr 



(3.18) 



(3.19) 



kh' Q - k 2 h x -™h + (l-l)(e + 2)Q = 2^GH 2 {rQ). 

(3.20) 

We eliminate then kho from Eqs. (3.19) and (3.20) to 
obtain 



-(1-H 2 r 2 ) 

\ ) dr 2 {l-H 2 r 2 ) rdr 



d 2 



(rQ) 



k 2 rQ 2d 

H -{rQ) 



-4ikH 2 rhi +2(1- H 2 r 2 ) H 2 r^- - 2k 2 K 

dr 



= -6H 2 (1 - H 2 r 2 ) h 



2ik (1 - H 2 r 2 ) - 2ikH 2 rh x - k 2 h 
v ' dr 



(3.28) 



+ 



(1 - H 2 r 2 ) 



2 r 



d ( 9 dh\ „ d ( 9 dK 
a ' 2— r 2 



dr \ dr J dr \ dr 



+ {1-H 2 r 2 ) 



£{e + i)h 



AH 2 r— + 2H 2 v 
dr 



dK_ 

dr 



+ {(, -!){(, + 2)^- = 2AirGH 2 {rQ). (3.21) 
r 

From Eq. (3.15) we finally obtain the wave equation 
d 2 Q 1 



where 

Veff 



t+l) 



(l - ffV) 



- 24irGH 2 



is an effective potential, and 

*° - (i) I <•«>■ 



(3.22) 



(3.23) 



(3.24) 



The derivation of the wave equation for the case of the 
even-parity perturbations has more algebraic complica- 
tions. From 5R^ V = Ah^ and the general expression 



= 6H 2 {1 - H 2 r 2 Y h 



(3.29) 



^{r 2 K)-2rh 
dr 



(1 - H 2 r 2 ) 



-i(l+l)K -2ikrh 1 = -6H 2 r 2 K, (3.30) 

where we have used Hi = Ho = h, which is allowed by 
the high symmetry of the de Sitter space, and Hi = hi. 
The three first-order equations (3.25)-(3.27) can be used 
to derive any of the subsequent second-order equations 
(3.28)-(3.30), provided the following algebraic relation- 
ship is satisfied 

[{£ -!){£ + 2)- AH 2 r 2 (l - 3H 2 r 2 )] h 



+ 



'+!) + 



2 {k 2 r 2 + 1 + 2H 2 r 2 ) 
(1 -H 2 r 2 ) 



K 



5 



2ikr + 



3iH 2 rl 



1) 



h! = 0. 



(3.31) 



The wave equation in a single unknown can now be 
derived from Eqs. (3.25), (3.26), (3.27) and (3.31). After 
a lot of algebraic manipulations we finally obtain 

cps i r 2 4(i-ffV) 

I AH r H - 

dr* 2 2 I r£>(r) 



dr* 



x [4H 2 r 2 (1 - 3H 2 r 2 ) -{I- + 2)] } 

if (1 - H 2 r 2 ) 

-Ak 2 + 4H 2 (2 + H 2 r 2 ) - A ^ - l)f + 2) 



4 4(1- H 2 r 2 ) 
_i_ + liL^fl i [ 4 ff 2 (1 - JfV) (1 - 3ff 2 r 2 ) + 2fc 2 



D(r) 
(1 - ffV) 



(£-!)(£ + 2) + 



3H 2 



1) 



5 = (l - H 



2 r 2 ) 



5 = 0, 
(3.32) 

(3.33) 



in which we have introduced the definition 

hi 
r 

and the function D(r) has the form 

2J.2 2 

D(r) = ( i-j?v) ~ 2(£ ~ 1)(f+2) 

+ (1 6g X 2 ) +4gV ( 1 ~ 3gV )- ( 3 - 34 ) 

We note that -D(r) approaches infinity as r tends to the 
cosmological horizon. The wave equation (3.32) can also 
be expressed in terms of another function S = h\/r 

d 2 S 2{l-H 2 r 2 ) 



dr* 



rD(r) 



x [4#V (1 - 3H 2 r 2 ) -(£- l)(t + 2)] 



dS 

dr* 



1 { (l-H 2 r 2 ) 
+ i 3F 2 (3 + H 2 r 2 ) - [ 2 > (t -!)(£ + 2) 



. 4 , 4(1 H 2 r 2 ) 



+k 2 -- + 



D(r) 



AH 2 (l-3H 2 r 2 ) 2 + 2k 2 



(Iz^_ 1)(<+2) + ^ (f+1) 



5 = 0, 



(3.35) 

with D{r) also as given by Eq. (3.34), which will be most 
useful in what follows. 



B. Odd-parity perturbations 

The stability of the static region of dc Sitter space to 
the odd-parity perturbations will be now examined. As 
in the Schwarzschild case [30,31], we shall distinguish two 
different cases: we first consider the situation when the 
frequency k is pure imaginary, and then we analyse the 
wave problem that results when that frequency is kept 
real. 

Set k = ia. Then, from Eqs. (3.22) and (3.24), we 
have 



= j(« 2 + Vi w )O = o 



(3.36) 



(3.37) 



The coordinate r* ranges from — oo to iir/H. The upper 
limit can of course be made zero by a re-definition of r* , 
such that now 



1 

H 



In 



1 - Hr 



1 + Hr 



The effective potential V e ff is real and positive every- 
where along this range, and vanishes at r* = — oo, i.e., 
on the boundary at r = H . The asymptotic solution 
to Eq. (3.36) as r approaches the cosmological horizon is 



Q;: ~ >' V P [ ± ^ ar " 



(3.38) 



As r — > 0, i.e., as one approaches the other boundary 
where we should require the perturbations to fall off to 
zero, one can write Eq. (3.36) in the form: 

^ + I(24^-ffi±il-o=) „,0, ,3.39, 

whose general solution can be given in terms of the Bessel 
functions [33] 

Qo - V~rC iVimjri Q VH 2 a 2 r) . (3.40) 

Since we should require the pertutbations to vanish as 
r — > 0, we have to choose for the Bessel function C = J 
[33]. However, if we take the function Q to be positive, 
one can see from Eq. (3.36) that d 2 Q/dr* 2 can never 
become negative within the entire range of r, from to 
H , and the solution that goes to zero at the origin r = 
(that is Eq. (3.40) with C = J) cannot be matched to 
the solution that goes to zero on the horizon at r = H^ 1 
(that is Eq. (3.38) with the sign + in the argument of 
the exponential.) It follows that the asymptotic solution 
near r = H^ 1 has to be 



(3.41) 
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where A is an arbitrary constant, so the radial solution 
near r = H~ 1 , Eq. (3.37), becomes 



h = - 



2H' 



(3.42) 



We can now compute the Kruskal perturbations given 
by Eqs. (3.13). Let us concentrate on the expression that 
results for h$\ near the surface r = H^ 1 (i.e., on u = v). 
Suppressing all angular dependence, we have 



[uh - v(l - H 2 r 2 )hi] e 



at 



H(u 2 -v 2 ) 



(^j+rv) Qh 



Aexp [-a(r* + t)] . 



Note that the sign of time has been reversed in the above 
equation with respect to those perturbations that appear 
in the Schwarzschild space. This reversal expresses the 
fact that the expansion of the universe is in many ways 
similar to the collapse of a star, except that the sense of 
time is reversed [21]. 

Using then exp(— r*) = l/{u 2 - v 2 ) 1 / 11 and e* = 

, _ ,1/2H 

I ^-j-ii J , and taking into account that u = v on the 
surface r — H -1 , we finally get 

&oa « \{-l) 1+a/2H A(u - v)-^+ 5a l 2H \u + v)- a ' 2H . 

(3.43) 

Now, since at t = (i.e. on u = 0) this perturbation 
becomes 

K^^{-l)- 2 ' H Av-^ H \ (3.44) 

which is clearly divergent at t = on the horizon, and 
the physically allowable perturbations should be regular 
everywhere in space at t = 0, we see that this perturba- 
tion is physically unacceptable, and hence cannot exist. 
It follows that the odd-parity perturbations with purely 
imaginary frequency ought to be ruled out. 

Let us consider now the solutions that correspond to 
real frequencies k. We shall look first at the case of in- 
going waves for which the asymptotic solutions near the 
horizon at r = H^ 1 are (V e ff = 0): 



Qh = Aexp 



ik 
2r* 



-oo. 



For these solutions it holds 

Aexp(-i/c/2r*) 



hi = 



H(l- H 2 r 2 ) 



(3.45) 



(3.46) 



ho = —cxp(-ik/2r*) . 



(3.47) 



For the cosmological perturbations in Kruskal coordi- 
nates we have then, 



ho3 oc 



A(u/2-v) 
H(u 2 - v 2 ) 



cxp (—ikr*/2 + ikt) . 



We note that for ingoing waves the horizon r 
should be taken at u = —v. Therefore, 



H- 



*03 



3^ 



4H(u 



■exp[ik(t-r*/2)}. 



Using again exp(— r*) = l/(u 2 — v 2 ) x l H and e* 
/ _ \^ 2H 

( ^Tv ) ' we nnau y obtain for this perturbation 



ho\ oc 



3A 



[i(u + v)]~ 



ik/H 



(3.48) 



AH(u + v) 

Since u + v — > on the horizon, at first glance this 
perturbation appears to be seriously divergent. However, 
as it also happens in the Schwarzschild case [30] , one can 
build wave packets which are convergent everywhere in 
space out of the monochromatic waves. If we form the 
purely ingoing waves into a wave packet by taking A to 
be a function of k given by the Fourier transform of a 
function f(q) = J A(k) exp(~ikq)dk, which vanishes for 
q < 1, by integrating over k, Eq. (3.48) transforms into 



ho\ a 



(3.49) 



There cannot be any singularity from the (u + v)^ 1 factor 
in Eq. (3.49) because / is nonzero only when i(u+v) > 1. 
Thus, /iQ3 does not diverge anywhere in space, but it is 
always pure imaginary. 

For outgoing waves, the asymptotic solutions near r = 
H^ 1 are given by (u = v): 

Q H = Acxj>{+ik/2r*) , r* -> -oo, (3.50) 

and, for wave packets formed as before, we finally obtain 

3 M 



h, 



AH(u-v) J [H 



f < — In [i(u 



(3.51) 



which, although coverging everywhere in space, is always 
pure imaginary, such as it happened for perturbation 
(3.49). The reason for these perturbations to be imag- 
inary resides in the fact that the argument of the log- 
arithm in the above expressions should be larger than 
unity, which in turn requires that both u and v are imag- 
inary simultaneously. Note that this does not imply rota- 
tion of the time t into the imaginary axis because Kruskal 
coordinates u and v appear in the form v/u in Eq. (3.12). 

Thus, the main conclusion for odd-parity perturba- 
tions is that, either they cannot exist whenever their fre- 
quency is purely imaginary, or they are stable and phys- 
ically acceptable if their frequency is real. Since in the 
latter case the solutions are purely imaginary, one can 
also conclude that such perturbations would not be ob- 
servable. 
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C. Even-parity perturbations 

In terms of the most covenient function S = 5/(1 — 
H 2 r 2 ), the second-order wave equation for the even- 
parity perturbations given by Eq. (3.35), near the hori- 
zon at r = if -1 , becomes 



d 2 S 1 



dr 



^ + -(fc 2 + 8ff 2 )S = 0. 



(3.52) 



Thus, the asymptotic form of the general solution at r = 
H~ x , Sh, would read: 



S H ~exp\ ±i\j-k 2 + 2H 2 r* 



(3.53) 



As for the boundary at r = 0, one attains that no per- 
turbation can consistently be expected, since from Eq. 
(3.35) we obtain So = (unless for t = for which case 
S is an arbitrary constant), and therefore h = hi = 0. 
Hence, we can readily get an expression for /iqo, indepen- 
dently of the value of the frequency. It follows that the 
two signs involved in the exponent of Eq. (3.53) arc al- 
lowed and should therefore be taken into account in our 
analysis. 

Again we first consider the case where the frequency k 
is purely imaginary, k = ia. Then 



Sjj oxp i 111- - ^a 2 r* 



(3.54) 



We have two distinct cases: case I, for which |a| > a c = 
2\/2H, and case II, for which \a\ < a c . In case I the 
asymptotic solution reduces to 



S H ~ exp (±£r*) , 



where 



2H 2 



(3.55) 



(3.56) 



is real. Because case II is qualitatively the same as that 
for real frequency (to be dealt with later on), we shall 
concentrate now on case I only. Since S — hi/r, near 
r = -ff -1 , one can assume the asymptotic forms 

hi = Acxp(±£r*) , h = 5exp(±£r*) . 

Choosing first the minus sign in the exponent of these two 
functions, from the equation relating radial functions h 
and hi to one another (which can be obtained by suitably 
combining Eqs. (3.26) and (3.27)), after specializing to 
r = H^ 1 and k — ia, i.e. 



AH 



dhi ,d 2 hi 



dr" 



+ 4 



dr* 



+ a 2 hi + 2a [ Hh + 2 



dh 
dr* 



B = 



2Hi + 4H 2 - a 2 
a(2£-H) 



A. 



(3.57) 



Let us now denote a — ±eH, with e > 2\/2- Then the 
analysis of all possible resulting cases (including the use 
of both signs in the exponent of Eq. (3.54)) will require 
considering 



2H. 



we get 



We shall look at three significant values of e, namely, 
2\/2, 3 and oo, in the following cases: (i) a > 0, £ > 
(B = V2A for e = 2y/2, A = for e = 3, and B = -A 
for e -> oo), (ii) a > 0, £ < (B = V2A for e = 2y/2, 
and A = B for e = 3 and £ — > oo), (iii) a < 0, £ < 
(B = -V2A for e = 2^2, B = -A for e = 3 and e — > oo), 
and (iv) a < 0, £ > (B = -V2A for e = 2\f2, A = 
for e = 3, and A = B for e — > oo). Clearly, if all of these 
particular cases led to stability of de Sitter space against 
the considered perturbations, one could conclude that de 
Sitter space is stable to such perturbations in all cases. 
Thus, for case I, the relation between the coefficients A 
and B for the asymptotic solutions h and hi must run 
between the extreme values for A (B fixed), A = B and 
A = —B, passing on A = 0. 

For A = B, the perturbation in the Kruskal coordi- 
nates, e.g. /jgo, is given by (angular dependence sup- 
pressed) 

h" F2{U ~ V) \ c-<* 

""00 — 9 9 " ,e 

u — v z 

= /n-a/a^^^LJy . (3.58) 

(u + w) 1+ V« 2 /4H 2 -2+ Q /2H 

At t = (u = 0), 

h 00 = (-l) 1 "\/ a2 /4H 2 -2-Q/H w -(2V" 2 /4H 2 -2+a/H) 

(3.59) 

Eq. (3.59) is divergent as v — > 0, except for the case 
a = — eH, £ > or £ very large, for which case /Iqo = 
—F 2 A. Hence, except for this case, all perturbations are 
physically unacceptable, as they all diverge at the initial 
time t = 0. 

For B = —A, we have 

I , \l-V" 2 /4H 2 -2+a/2ff 

Co = (-l)-^ 2H F 2 A {U + V) V 7 , 

( u _ w )l + V« 2 /4H 2 -2+a/2H 

(3.60) 

which at t = (m = 0) reduces to 

(3.61) 
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We note that all of these perturbations are physically 
unacceptable, except for the case a = — eH, £ < 0, with 
e = 3, where = iF 2 Av. 

Finally, when A — 0, we obtain 



h. 



K 
00 



= {-l)- a ' 2H F 2 B(u 2 + v 2 ) 



(u + w )«/2ff-l-V« 2 /4ff 2 -2 
( u _ v y + ^/ a y4H2-2+a/2H ' 

(3.62) 



Again at t = (u = 0), this perturbation reduces to 

fc* = (_l)-( 1 +V« 2 /4H 2 -2+a/H) F 2^-2Va 2 /4^-2_ 



It can easily be checked that in all the cases, without any 
exception, Eq. (3.63) diverges as v — ► 0, and therefore 
this perturbation is physically unacceptable and should 
be ruled out. 

We are in this way left with two physically acceptable 
even-parity perturbations for purely imaginary frequency 
in case I: that given by Eq. (3.58) for negative a, positive 
£ and very large e, and that given by Eq. (3.60) for neg- 
ative a, negative £ and e about 3. Since negatives values 
of a correspond to the case of outgoing perturbations for 
which u = v on the horizon r = H -1 , these perturba- 
tions will be stable as the resulting powers to the factor 
(u — v) are positive definite in both cases. Note, further- 
more, that at least the perturbation given by Eq. (3.60) 
is always purely imaginary. Thus, even-parity perturba- 
tions with purely imaginary frequency in case I either are 
physically unacceptable or, unlike in the Schwarzschild 
space [30] , are stable and most of them imaginary. 

In the case that k is kept real the asymptotic general 
solution on the cosmological horizon has already been 
given by Eq. (3.53). The analysis to follow will also 
be valid for purely imaginary frequency satisfying the 
condition implied by case II, with the asymptotic solution 
being given by Eq. (3.54) in this case. For real frequency, 
the relation between h and hi is 



k 2 h x -A ff— i 



jjdh\ d 2 hi 
dr* dr* 2 



2k 



„ dh 



(3.64) 



where we have specialized to the region r = H . Let 
us now assume, like it was made for purely imaginary 
frequencies, that the asymptotic forms for h and hi near 
the cosmological horizon are given by 

hi = Aexp(— iur*), h = i?exp(— ivr*), (3.65) 

in which we have introduced the short-hand notation 



v = \j2H 2 + -k 2 . 



(3.66) 



From Eqs. (3.64) and (3.65) we get 



B 



+ 



4iH 3 



2v (k 2 + 5H 2 

k{k 2 + 9H 2 ) ' k(k 2 + 9H 2 ) 



A. (3.67) 



Restricting to the case k 2 >> H 2 , so that v ~ k/2, we 
see that there are two solutions: when k ~ 2v (ingoing 
waves), A = —B, and when k ~ — 2v (outgoing waves), 
A = B. Suppressing again the angular dependence of the 
perturbations, we can then compute such perturbations 
in Kruskal coordinates. In the case of ingoing waves, 
u = — v on the horizon, and we have 



h^ = F 2 A(-iy k ' 2H 



(u + v) 1 



ik/H 



u — V 



(3.63) and for outgoing waves (u 



h& = F"A{-l) 



ik/2H 



(u — v) 



1+ik/H 



U + V 



(3.68) 



(3.69) 



Clearly, the perturbations (3.68) and (3.69) are stable 
near the cosmological horizon and everywhere inside it, 
even in the forms given by these equations, without build- 
ing suitable wave packets by superposing them. This 
analysis can readily be generalized to any values of k 
and H, obtaining the same conclusion. 



IV. STABILITY OF MULTIPLY CONNECTED DE 
SITTER SPACE 

In Sec. Ill we have investigated the stability prop- 
erties of the simply connected de Sitter region covered 
by static coordinates. We had to do so because, as far 
as I can know, such an analysis had not been carried 
out so far, and we needed it to prepare our system to 
study the perturbations when the coordinates involved 
are identified in such a way that this region of de Sitter 
space becomes multiply connected topologically, such as 
it was discussed in Sec. II. We have obtained that on the 
static region, the de Sitter is also stable to the first-order 
perturbations that satisfy its symmetries. 



A. Classical perturbations 

We shall now study the effect that topological multi- 
ple connectedness has on the stability of de Sitter space. 
Because of the high symmetry of this space, the time pa- 
rameter t always appears in the form of a factorized expo- 
nential factor in all the perturbations, cither as exp(±crf), 
if the frequency of the perturbative modes is pure imagi- 
nary, or as exp(±ifct), if that frequency is kept real. Thus, 
we can generically write the time factor as exp[<7(fc)i], 
with g{k) a given function of the mode frequency. From 
our discussion in Sec. II, it follows that de Sitter space 
can be made multiply connected by simply including the 
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time identification t <-> t + nb/H, with 6 a dimcnsionless 
arbitrary period and n any integer number. This will 
amount to the insertion of an additional factor 

exp [nbg(k)/H] 

for each value of the integer n in the distinct expressions 
for the first-order perturbations obtained in Sec. III. In 
order to take into account all possible values of n, along 
its infinite range, one then should coveniently sum over 
all n, from to oo, introducing the statistical factor 1/n! 
to account for the equivalent statistical weight one must 
attribute to all of these contributions. Thus, the general 
expression for the perturbations in multiply connected de 
Sitter space would be 



h?Ab) 



E 

n=0 



h$exp [nbg(k)/H] 



(4.1) 



where hfj generically denotes the first-order perturba- 
tions in Kruskal coordinates for simply connected de Sit- 
ter space which were computed in Sec. III. 

In what follows we shall perform the calculation of 
the relevant hfj(b) for all physically acceptable perturba- 
tions. We shall first restrict ourselves to the regime where 
both kb/ H and ab/H are much smaller than unity; i.e. 
we will work in the regime characterized by nonchronal 
regions and CTC's whose size is very small. As it will 
be seen below, this is the regime of most physical in- 
terest where vacuum quantum fluctuations can be kept 
convergent everywhere. Let us start with odd-parity per- 
turbations with real frequency. For the case of ingoing 
waves, we have for the asymptotic solution at r = H^ 1 
[34] 



"03 



(b) = 



3A 



AH(u + v) 



[i(u + v)]' 



00 Ankb/H 
ik/H ^ e 

n=0 



ZA 



AH{u + v) 



exp 



-ik\n[i(u+v)]/H 



(4.2) 



For small values of kb/H, Eq. (4.2) can be approximated 
to: 



3Ae 



AH{u + v] 



■ exp 



H 



[\n(i(u + v))-b}\ . 



(4.3) 



Forming again a wave packet out of monochromatic per- 
turbations (4.3), we finally obtain for this type of waves 



3e 



AH{u + v) 



f 



1 hn(u + v) + y - b 



(4.4) 



which still is a convergent expression for all times. If we 
let b to be complex, so that b = b + in/2, then Eq. (4.4) 
becomes not only convergent but pure real as well. 



For outgoing waves, an analogous calculation leads fi- 
nally to the perturbation: 



3e 



AH(u - v) 



f 



1 hn(u - v) + y + b 



(4.5) 



which is also always covergent and pure real if we sim- 
ilarly let b be complex and given by b = b — i-k/2. We 
can then conclude that in the considered regime, multi- 
ply connected de Sitter space is stable to all odd-parity 
perturbations that are physically acceptable. 

For even-parity perturbations which also are physically 
acceptable we have [34] 



-nba/ H 



n=0 



ni 



- h K 



oo ex P 



-ba/H 



(4.6) 



for any value of kb/H. Thus, making the de Sitter space 
multiply connected preserves the stability of these per- 
turbations and increases their amplitude, specially for 
small values of ba/H. 

Next we consider even-parity perturbations with real 
frequency. We first note that in this case the pertur- 
bations corresponding to the asymptotic solutions near 
r = H^ 1 can also be expressed in terms of wave packets 
in the simply connected case. They are: 



Co = 



u — v 



f 



for ingoing waves, and 



u + v 



1 f, / \ in 
— \ln(u + v) + — 



-hn(u-v) + - 



(4.7) 



(4.8) 



for outgoing waves. Because of the form of the Kruskal- 
coordinate dependent prefactor, these expressions arc 
real in any case. 

When we multiply connect the de Sitter space in the 
regime of small values of kb/H, Eqs. (4.7) and (4.8) 
transform into: 



h«{b) - eF 2 



u + v 



f 



1 

H 



ln(u + v) + 



in 



(4.9) 



for ingoing waves, and 



h&(b) = eF 2 



u + v 



— [ ln(u 



(4.10) 



for outgoing waves. Note that the argument of the func- 
tion / becomes real when we allow b to be complex and 
given by b — b ± iir/2, with the + sign for ingoing waves 
and the - sign for outgoing waves. Anyway, the per- 
turbations given by Eqs. (4.9) and (4.10) keep being 
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convergent and, therefore, one can conclude that in the 
regime of very small values of b/H, the multiply con- 
nected de Sitter space is also stable to first-order pertur- 
bations which respect the symmetry of this space. 

As to the perturbations for larger values of b/H and 
real frequency, we first note that the expressions for the 
components fiQ 3 (b) and hff (b), before forming any wave 
packets, are given by expressions which are the same as 
those obtained above in the regime of very small values 
of b/H, but with the parameter b replaced for the k- 
dependent function —^jfexp^kb/H). The wave packets 
formed using the same procedure as for all the above 
cases will therefore involve characteristic functions of the 
form 



'{!/ 



dxbe x 



b + In [i(u ± 



«)]}' 



where x = exp(ikb/H), and the sign +/- in the argument 
of the logarithm stands for ingoing/ outgoing waves. Now, 
for b >> In \i(u ± v)\ we can obtain [33,34] 



H 



Ei{x)\l 



with Ei{x) the exponential integral function. Since the 
argument of / is then always smaller than unity, we have 
/ = on this regime. 

As (u ± v) becomes very small, so that b « y = 
|Reln [i(u ± v)] |, the wave-packet function approaches 
the form [33,34] 



y^oo 



+ 1;1) 



y 



o, 



with $ the degenerate (confluent) hypergeometric func- 
tion. Thus, also for arbitrarily large nonchronal regions, 
the multiply connected de Sitter space is stable to all 
physically allowable classical perturbations. 



B. Quantum fluctuations 

In what follows we shall briefly discuss the possible 
influence that multiply connectedness may have on the 
quantum stability of the de Sitter space. Because of the 
presence of a chronology horizon on the surface r = H -1 , 
it could at first sight be thought that the quantum renor- 
malized stress-energy tensor, (T AI ^) ren , for vacuum quan- 
tum fluctuations generated in multiply connected de Sit- 
ter space ought to diverge [35]. However, it has recently 
been stressed that this could not be actually the case if 
either we consistently redefine the quantum vacuum [23] , 
or we introduce a suitable quantization of the relevant 
time parameter, beyond semiclassical approximation [27]. 
To see how these ideas apply to the case under study, let 
us work in the Euclidean framework where the Kruskal 
metric is obtained by rotating coordinates u and t to the 



imaginary axis, starting with Eq. (3.9); i.e., u — ir] and 
t = it. We get then 



ds 2 = 



H 2 (l + T] 2 +V 2 



(drj 2 + dv 2 ) +r 2 dfll, (4.11) 



which in fact is definite positive. Corresponding to this 
metric, Euclidean time r will be defined by the relation 



cxp(2iHr) 



IT] 



(4.12) 



Wick rotating also in the identification t <-> n/3 (where the 
nonperiodic time term is disregarded and (3 — b/H is the 
period) that makes de Sitter space multiply connected, 
Eq. (4.12) transforms into 



exp(2inH/3) = 



V — IT) 



V + 17] 

from which one can get the complex relation 
v — irj = \Jv 2 + r] 2 exp(inHf3). 



(4.13) 



(4.14) 



It follows from Eq. (4.14) that the Euclidean time pre- 
serves a periodic character also on the Euclidean sector: 
(3 = 2n/H, that is to say, b = 2ir. This result can be 
interpreted along the following lines. First of all, one can 
readily see that multiply connectedness in de Sitter space 
is nothing but the Lorentzian counterpart of the existing 
thermal states that are uncovered in the Euclidean de- 
scription [36] . This relation might be reflecting the origin 
of the excess of some perturbative waves which we have 
found above for multiply connected de Sitter space with 
respect to its simply connected space. 

Quite more importantly, the value b — 2n can be used 
(as Gott and Li did [23]) to redefine a conformal vacuum 
in Euclidean space for which (T Miy ) re „ does not diverge 
even on the chronology horizon. However, the meaning 
of this horizon in such a vacuum has been discussed by 
Kay, Radzikowski and Wald [24] and Cassidy [25], so 
that some quite founded doubts can be cast on its real 
existence and capability to restore quantum stability this 
way. But if we adhere to the also recently suggested [27] 
kind of time quantization by which t = (n+j)to (with to a 
constant time whose value is of the order the Planck time, 
and 7 the automorphic parameter [37,38], < 7 < 1/2), 
and note the formal analogy of this expression with that 
is implied by the identification t <-» t + nb/H when we 
take b — 2n and t = 2-kj/H, we see that quantum sta- 
bility could be unambiguously restored in multiply con- 
nected de Sitter space, provided we accept restricting the 
nonchronal region and the CTC's on it to be essentially 
at the Planck scale, 1/H - l p [27]. 



V. SUMMARY AND CONCLUSIONS 

The main aim of this paper is to study the classical 
and quantum stability of the multiply connected de Sit- 
ter space. This space arises when we introduce some 
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periodicity conditions on the coordinates describing the 
five-dimensional Minkowski hyperboloid and can only be 
exhibited on the region covered by the static de Sit- 
ter coordinates. Using a first-order perturbation formal- 
ism analogous to that which was originally developed by 
Regge and Wheeler for Schwarzschild metric [26] , we have 
shown that multiply connected de Sitter space is classi- 
cally stable to these perturbations, no matter the size of 
the static region. Although stability aginst higher-order 
perturbations has not been checked in this paper, one 
would expect these perturbations not to introduce any 
instabilities, such as it occurs in Schwarzschild spacetime 

By continuing the Kruskal extension of the multiply 
connected, static de Sitter metric into its Euclidean sec- 
tion, we have also argued that quantum vacuum fluctu- 
ations should not induce any divergence on this space, 
provided the nonchronal region and the CTC's on it are 
all sufficiently small, probably of the order the Planck 
size [27]. We therefore consider multiply connected de 
Sitter universes to be genuine components of any future 
description of a well-defined theory of quantum gravity. 
In particular, the considered stable little multiply con- 
nected universes should be included, together with Eu- 
clidean and multiply connected wormholes [19], ringholes 
[39], Klein bottleholes [40] and virtual black holes [41], 
as components of the vacuum quantum spacetime foam, 
where their CTC's would contribute the required viola- 
tion of causality that governs the foam. Thus, Planck- 
sized de Sitter universes containing CTC's can help to 
define the boundary conditions of the universe we live in, 
probably along the lines recently suggested by Gott and 
Li [18] and discussed in the Introduction of this paper. 

To close up, I would like to refer to the interesting pos- 
sibility that the Hartle-Hawking and the Gott-li condi- 
tions might both be seen to imply initial physical pictures 
which, at least in a way, appear to be complementary. In 
classical relativity time and spatial coordinates can still 
be distinguished by the fact that, whereas spatial coor- 
dinates do not single out a particular sign to run over, 
time can only run forwards. This does not help, none the 
less, to understand the Lorentzian signature of the clas- 
sical metrics which are not positive definite. However, 
as we approach the regime of the quantum spacetime 
foam, all of this distinction can be thought to vanish, 
since in such a regime there will be CTC's everywhere 
and hence the two time direction would become equally 
allowable and, at the same time, metrics can be taken 
to be positive definite. In order to describe the quantum 
origin of the universe, one then can either keep CTC's 
and Lorentzian signatures simultaneously, as assumed by 
Gott and Li, or disregard CTC's while using Euclidean 
signature where time becomes spacelike, as suggested by 
Hartle and Hawking. As seen in this way, the two pictures 
would actually describe rather complementary aspects of 
the initial physical situation. 
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